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Abstract 

An algorithm recently presented by Lake to obtain all static spher- 
ically symmetric perfect fluid solutions, is extended to the case of lo- 
cally anisotropic fluids (principal stresses unequal). As expected, the 
new formalism requires the knowledge of two functions (instead of 
one) to generate all possible solutions. To illustrate the method some 
known cases are recovered. 
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1 Introduction 



As is well known, static spherically symmetric perfect fluid distributions in 
general relativity, are described by a system of three independent Einstein 
equations for four variables (two metric functions, the energy density and the 
isotropic pressure). Thus, aditional information in the form of an equation of 
state or an heuristic assumption involving metric and/or physical variables 
has to be provided in order to integrate the system. This situation suggests 
the possibility of obtaining any possible solution, giving a single generating 
function. A formalism to obtain solutions in this way has been recently 
presented by Lake [1] (see also [2]). 

The purpose of this work is to extend the above mentioned formalism to 
the case of locally anisotropic fluids. 

The motivation for doing so is provided by the fact that the assumption of 
local anisotropy of pressure, which seems to be very reasonable for describing 
the matter distribution under a variety of circumstances, has been proved to 
be very useful in the study of relativistic compact objects (see [3]- [13] and 
references therein). 

In the next section we shall present the general equations and the formal- 
ism to obtain the solutions, then we shall apply the method to analyze some 
specific cases. 



2 The Einstein equations for static locally anisotropic 
fluids 



which has to satisfy the Einstein equations. For a locally anisotropic fluid 
they are 



In curvature coordinates the line element reads 



ds 2 = - 



e "(r) dt 2 + e X(r) dr 2 + ^2 + ^2^2 



(1) 






(3) 
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8nP ± = i^v" + v 12 - AV + , (4) 

where primes denote derivative with respect to r, and p, P r and Pj_ are the 
proper energy density, radial pressure and tangential pressure respectively 

2.1 The algorithm 

From (3) and (4) it follows: 

&r(P r - P,) = - (-)» + Tr + - 2 ) + e~>-(- + -)- 7r (5) 

Then, introducing the variables 

e "( r ) = e J (2z(r)-2/r)dr ^ 

and 

e" A = y(r) (7) 

and feeding back into (5) we get: 

r 2z' 6 4. 21 

y' + y — + 2z - - + — = — - + n r , 8 

z r r 2 z z r z 

with n(r) = 8n(P r -P ± ). 
Integrating (8) we obtain for A: 

e x(r) = . (9) 

r 6(-2 / z(r)(l+n(r)r 2 )e J rS I rfr + 

where C is a constant of integration. Then, using (6) and (9) in (1) we get: 

ds 2 = - e n 2z ^- 2 ^ dr dt 2 + ^ (r)e /(^+^(^ dr 2 +r 2 de 2 +r 2 sin 2 ed(t ) 2 . 

(10) 
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Thus any solution describing a static anisotropic fluid distribution is fully 
determined by means of the two generating functions II and z. 

It will be convenient to express the physical variables in terms of metric 
and generating functions, in order to impose conditions leading to physically 
meaningful solutions. Thus we have: 

AnP r = z(r-2m)+m/r-l (n) 



4n P = ^ (12) 



and 



2m. , , 9 z 1 . .mm!. 
4nP ± = 1 )( z ' + z 2-- + -)+z(- 13 

where the mass function m(r) is defined as usual by 



g _ A = 1 _ Mr) 
r 

Physically meaningful solutions must be regular at the origin, and should 
satisfy the conditions p > 0, p > P r ,P±. If stability is required then p and 
P r must be monotonically decreasing functions of r. 

To avoid singular behaviour of physical variables on the boundary of 
the source (S), solutions should also satisfy the Darmois conditions on the 
boundary . Implying (P r )s = and 

= e~ As = 1 - — (15) 

with m 2 = M, and r 2 denotes the radius of the fluid distribution. 

2.2 The locally isotropic case 

If we impose the isotropic condition on pressure 

n = 8vr(P r - P ± ) = (16) 

in (10) we obatin: 
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ds 2 = - e f(Mr)-2/r)dr dt 2 + g W dr 2 +rW+r W## 2 

r 6(_ 2 ; rz( 8 r rf r + C) 

(17) 

which is the same result obtained in [l],with z(r) = $(r)' + K 

3 Some examples 

We shall next apply the algorithm to reproduce some known situations. 
3.1 Conformally flat anisotropic fluids 

Instead of giving two generating functions, we may provide one generating 
function and an additional ansatz. Thus for example, in the spherically 
symmetric case we know that there is only one independent component of 
the Weyl tensor. Therefore the conformally flat condition reduces to a single 
equation which reads 

v" V. 2 v'\' v 1 — A' l-e A ft 

Y + <2> - — - — + — = ' (18) 

Equation (18) has been integrated in [14], giving: 

/C 2 
— dr), (19) 

which, in term of z becomes 

2 



+ —tanh([—dr). (20) 
r J r 

On the other hand from (4) and (18), it follows: 



z — — 
r 



n = r(i-A'. (21) 



Thus the system is completely determined (in this case) provided a single 
generating function z is known. 
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3.2 Bowers— Liang solution 



This solution corresponds to an anisotropic fluid with an homogeneous energy 
density distribution p = p = constant [15], and is given by: 

"3 (1 - 2M/r s f /2 - (1 - 2m/r) h/21 



/ \ 4vr 3 
m(r) = —r p 



Air 

M = —rlp 



The two generating functions for this metric are: 



z = 



r z v r ' I _ 

(1-^)4 r 



3(1-^)4 



and 



n 



-6C 



^_ 1)2 (1 _ 2M)# 



(22) 
(23) 

(24) 
(25) 



with /i = 1 — 2C = Constant. The case /i = 1 reproduces the well known 
Schwarzschild interior solution, whereas the case h = describes the Florides 
solution [16]. 



3.3 Anisotropic solutions with a non-local equation of 
state 

An interesting family of solutions may be found from the assumption that the 
energy density and the radial pressure are related by a non-local equation of 
state of the form [17] 

P r {r) = p{r) - - f r f 2 p(f)df + ^—r (26) 
r A Jo Airr 6 

or, using(12) 

, , ml m C 
r[r) = W " 2^ + 2^' { } 

From (11) and (27) it follows that these solutions are defined by the gener- 
ating function z of the form: 

rml — 3m + 2C + r . . 

z = 7 — ; • (28) 

r(r - 2m) v ' 
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